TO APPEAR IN The Astrophysical Journal. 

Preprint typeset using I^TfiX style emulateapj v. 26/01/00 



BULGE AND HALO KINEMATICS ACROSS THE HUBBLE SEQUENCE 
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ABSTRACT 

The correlation between the maximum rotational velocity of the disk (v m ) and the central stellar velocity dis- 
persion of the bulge (cto) offers insights into the relationship between the halo and the bulge. We have assembled 
integrated H I line widths and central stellar velocity dispersions to study the v m —ao relation for 792 galaxies span- 
ning a broad range of Hubble types. Contrary to earlier studies based on much smaller samples, we find that the 
v m -ao relation exhibits significant intrinsic scatter and that its zeropoint varies systematically with galaxy mor- 
phology, bulge-to-disk ratio, and light concentration, as expected from basic dynamical considerations. Nucleated 
but bulgeless late-type spiral galaxies depart significantly from the v m -ao relation. While these results render 
questionable any attempt to supplant the bulge with the halo as the fundamental determinant of the central black 
hole mass in galaxies, the observed distribution of v m /ao, which depends on both the density profile and kinematic 
structure of the galaxy, offers a useful constraint on galaxy formation models. With the aid of a near-infrared 
Tully-Fisher relation, we identify a population of otherwise normal, luminous galaxies that have exceptionally 
low values of v m /ao. We argue that a significant fraction of the H I gas in these kinematically anomalous objects is 
dynamically unrelaxed, having been acquired externally either through capture from tidal interactions or through 
cold accretion from the intergalactic medium. 

Subject headings: galaxies: bulges — galaxies: ISM — galaxies: kinematics and dynamics — galaxies: nuclei 



1. INTRODUCTION 

It is accepted that every galaxy contains an extended dark 
matter halo, but precisely how the halo couples to the lumi- 
nous components of the galaxy remains a subject of lively de- 
bate. Much attention has been given to the connection between 
the disk and the halo, in particular as it concerns the fraction 
of baryons that collapse to form the disk and the role of adia- 
batic contraction (e.g., Navarro & Steinmetz, 2000; Dutton et 
al. 2007). In disk galaxies that contain bulges, what is the re- 
lationship between the bulge and the halo? While a variety of 
methods can be used to address this question, important insights 
can be gained by investigating the relation between the line-of- 
sight central stellar velocity dispersion of the bulge, ctq, and the 
deprojected maximum rotation velocity of the disk, v m , which 
effectively traces the circular velocity of the halo. Whitmore 
and collaborators (Whitmore et al. 1979; Whitmore & Kirsh- 
ner 1981) first looked into this issue, using a small sample of 
spiral galaxies for which they had ern and v m measured through 
integrated H I velocity profiles. They find v m /ao as 1.2 to 2, 
with the ratio increasing with decreasing bulge-to-disk ratio; 
elliptical galaxies roughly follow the same pattern (Fall 1987; 
Franx 1993). These trends were subsequently confirmed for a 
larger sample of normal (Whittle 1992a) and Seyfert (Nelson & 
Whittle 1996) galaxies. 

The v m -<ro relation recently has attracted renewed attention 
in the context of black hole demographics studies. The interest 
is three-fold. First, if crn is related to v m , then the existence of 
the Mbh - o~q relation (Gebhardt et al. 2000; Ferrarese & Mer- 
ritt 2000) suggests that the black hole mass may be more funda- 
mentally tied to the halo mass rather than the bulge mass. This 
is the argument made by Ferrarese (2002), who revisited the 
v m -ao correlation originally introduced by Whitmore and col- 
laborators. Second, apart from its theoretical implications, the 
existence of a v m — <7n correlation presents a new empirical tool 



for black hole demographics studies. Stellar velocity disper- 
sions are not always easy or even possible to obtain, especially 
for active galaxies, whose bright nonstellar nuclei often over- 
whelm the stellar continuum and make measurement of cen- 
tral stellar velocity dispersions exceedingly challenging (e.g., 
Greene & Ho 2006). In such circumstances, it may be more 
feasible to measure v m for the disk, either through spatially re- 
solved rotation curves or integrated H I line profiles (Ho et al. 
2007a). Finally, the v m -<ro correlation represents a new scaling 
relation for galaxies, which, like other more familiar scaling 
relations, serves as an important boundary condition for theo- 
retical models of galaxy formation. 

Ferrarese (2002) compiled kinematic data for 16 disk (spi- 
ral and SO) galaxies with dynamical determinations of black 
hole masses to show that the rotation velocity of the disk, mea- 
sured on the flat part of the rotation curve, follows a tight, nearly 
linear correlation with the central stellar velocity dispersion of 
the bulge. For her sample, the correlation breaks down for cto 
<J 80 km s . The 19 elliptical galaxies with rotation veloci- 
ties derived from the dynamical models of Kronawitter et al. 
(2000) and Gerhard et al. (2001) seem to fall on the same 
correlation, suggesting that the v m -ao relation is universally 
obeyed by galaxies of all types. This result was echoed in sub- 
sequent studies (Baes et al. 2003), some of which (Pizzella 
et al. 2005; Buyle et al. 2006) additionally suggested that low- 
surface brightness galaxies follow a separate, nearly parallel re- 
lation compared to high-surface brightness galaxies by having 
a larger v m for a given o-q. In the most comprehensive anal- 
ysis to date, Courteau et al. (2007) challenged the existence 
of a tight, Hubble type-invariant v,„-o-q relation. Instead, their 
enlarged sample clearly demonstrates that v m /ao systematically 
varies with the concentration of the galaxy light profile, a result 
already apparent in the original work of Whitmore and others. 

This paper examines the v m —(7Q correlation using an exten- 
sive sample of nearby galaxies with accurate measurements of 
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TABLE 1: The Sample 
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0.44 
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68.2 


437 


4.5 


-19.11 


-23.28 


9.46 


ESO 186-055 


2.0 


N 




0.43 


1 


71.0 


4708 


64.5 


-20.02 


-23.04 




ESO 189-007 


3.8 


Y 




0.72 


1 


45.2 


2982 


39.5 


-22.31 


-24.03 


10.16 


ESO 206-014 


5.0 


Y 




0.83 


1 


34.5 


15000 


210.5 


-22.47 


-23.65 




ESO 215-039 


4.8 


Y 




0.71 


1 


46.2 


4348 


59.1 


-21.32 


-23.29 


9.66 



NOTE. — Col. (1) Galaxy name. Col. (2) Morphological type index. Col. (3) Presence of a bar. Col. (4) The concentration 
index in the i band, defined as C59 = 5 log(rpgo/rp5o), where rpgo and rpso are the Petrosian radii enclosing 90% and 50% of 
the light, respectively. Col. (5) Ratio of semi-minor to semi-major axis isophotal diameter measured at a surface brightness of 
fiB = 25 mag arcsec -2 . Col. (6) Reference for a and b. Col. (7) Inclination angle, calculated from b/a, as described in the 
text. Col. (8) Heliocentric velocity. Col. (9) Luminosity distance, derived assuming a Local Group infall velocity of 208 km 
s -1 toward the Virgo cluster and a Hubble constant of Ho = 70 km s _1 Mpc -1 . Additional data come from Rizzi et al. 2006 
(NGC 300), Lee et al. 2002 (NGC 598), Thim et al. 2004 (NGC 4395), and Karachentsev et al. 2003 (NGC 7793). Col. (10) 
Total absolute £?-band magnitude, corrected for Galactic and internal extinction, using the inclination angle from Col. (7) and 
the formalism of Bottinelli et al. 1995. Col. (11) Total absolute iG-band magnitude from 2MASS; these values have not been 
corrected for Galactic or internal extinction. Col. (12) H I mass, calculated using the following prescription from the RC3, 
log(M H i/M ) = -0.4m 2 i + log(l + z) + 21og(£>i/Mpc) + 12.3364, where m 2 i is the H I "flux" and 2 is the redshift. Data for 
Cols. (2), (3), (8), (10), and (12) are taken from HyperLeda. Table 1 is presented in its entirety in the electronic edition of the 
Astrophysical Journal. A portion is shown here for guidance regarding its form and content. 

REFERENCES.— (1) RC3; (2) Nilsen 1973; (3) Paturel et al. 2000; (4) Freudling et al. 1992; (5) Vorontsov-Velyaminov et al. 
1963; (6) MacGillivray et al. 1988; (7) Haynes et al. 1997; (8) HyperLeda; (9) Binggeli et al. 1985. 



both central stellar velocity dispersions and disk rotation ve- 
locities derived from integrated H I line widths. Our sample 
is significantly larger than those employed in recent studies, 
and it covers galaxies spanning a very wide range in Hubble 
types, allowing us to investigate trends with bulge-to-disk ra- 
tio. This analysis is similar in spirit to that of Courteau et al. 
(2007), the main difference being that ours utilizes integrated 
H I line widths instead of resolved rotation curves to estimate 
v m , a shortcut that enables us to dramatically increase the sam- 
ple size. 

We confirm, with greater statistical weight, that the v m -ao 
relation systematically depends on the galaxy luminosity den- 
sity profile, parameterized by Hubble type, bulge-to-disk ratio, 
or concentration index, and that there exists significant intrin- 
sic scatter at a given density. This finding does not bode well 
for attempts to use galaxy rotation velocities to predict black 
hole masses, nor does it support the contention that the black 
hole mass is more fundamentally linked to the halo than the 
bulge, but it does offer a new, potentially powerful constraint 
on galaxy formation models. We draw attention to a subset of 
galaxies characterized by having exceptionally narrow H I pro- 
files for their central stellar velocity dispersion. We suggest that 
the neutral hydrogen in these systems lies in a plane offset from 
that of the stellar disk, or that it has not yet settled into dynam- 
ical equilibrium with the stars. 

2. DATA COMPILATION 

Our goal is to compile as large a collection of galaxies as pos- 
sible having relatively homogeneous, well-documented mea- 
surements of central stellar velocity dispersions and rotational 
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velocities that probe the flat part of the disk rotation curve. 
Apart from the sample size, a key difference between our anal- 
ysis and those of previous studies [with the exception of ini- 
tial work by Whitmore et al. (1979), Whitmore & Kirshner 
(1981), Whittle (1992a), and Nelson & Whittle (1996)] is that 
our values of v m come from spatially integrated (single-dish) 
H I profiles. Previous authors have stressed the importance of 
measuring v m from spatially resolved, extended optical spectra 
that sample the flat part of the rotation curve. Here we wish 
to emphasize that integrated H I line widths provide a robust 
and efficient substitute for estimating v m . This has been well- 
documented in numerous studies, many motivated by the desire 
to use H I line widths for distance-scale investigations using the 
Tully-Fisher (Tully & Fisher 1977) relation (e.g., Rubin et al. 
1978; Thonnard 1983; Mathewson et al. 1992; Courteau 1997). 
The global H I line width imprints both the shape of the galaxy's 
rotation curve and the actual spatial distribution of the neutral 
hydrogen. But since the H I distribution in spiral galaxies typi- 
cally extends to twice the optical radius (Broeils & Rhee 1997; 
Noordermeer et al. 2005), in practice the width of the H I ve- 
locity profile is quite robust to different rotation curves and H I 
distributions (Roberts 1978). 

We have compiled a database that consists of four samples, 
which we describe in turn. 

• Sample 1 — The primary source of data for our inves- 
tigation comes from Hyperleda^ (Paturel et al. 2003a). 
This catalog is periodically updated, and the entries used 
in the present analysis are reported to be current up to the 
end of 2003. The stellar velocity dispersions in 
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TABLE 2: Kinematical Data 
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ESO 097-013 
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18.8 
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115.1 


8.8 


14.5 


10 
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2.0 
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11.0 










ESO 189-007 


4a 


91.0 


2.0 


317.0 


295.1 


178.7 


14.7 


14.5 


10 


1 




ESO 206-014 


4c 


54.0 


2.0 






199.9 


5.0 










ESO 215-039 


4a 


71.0 


11.0 


267.9 


215.8 


132.2 


4.7 


14.5 


10 


2 





NOTE. — Col. (1) Galaxy name. Col. (2) Sample. Col. (3) Central stellar velocity dispersion. Col. (4) Error on cq. Col. (5) 
Width of the H I profile at 20% of the maximum peak. Col. (6) Width of the H I profile at 50% of the maximum peak. Col. (7) 
Maximum rotation velocity v m , corrected for instrumental resolution, interstellar turbulence, (1 + z) cosmological stretching, and 
inclination angle. We obtained v m from either W20 or ITso using the calibration given in Paturel et al. 2003b; when both IT20 
and W50 are available, v m is the average of the two derived values. Col. (8) Error on v m , assumed to be approximated by the 
error on the homogenized value of the rotation velocity given in HyperLeda. Col. (9) Radius searched for interlopers and confusing 
sources. Col. (10) Telescope used to obtain the H I data: (1) Nangay; (2) Westerbork Sythesis Radio Telescope; (3) NRAO Green 
Bank 140 foot: (5) Jodrell Bank 76 m; (7) NRAO Green Bank 300 foot; (8) Effelsberg 100 m; (9) Arecibo; (10) Parkes 64 m: (14) 
Jodrell Bank 125 footx85 foot; (15) Very Large Array. Col. (11) Reference for the H I data. Col. (12) Notes: (1) Faint galaxy in 
the beam, but unlikely to be a source of confusion. (2) Significant galaxy in the beam, but not a source of confusion according to 
interferometric data. Table 2 is presented in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown 
here for guidance regarding its form and content. 

REFERENCES.— (1) Koribalski et al. 2004; (2) Mathewson & Ford 1996; (3) Paturel et al. 2003b; (4) Theureau et al. 1998; (5) 
Mould et al. 1991; (6) Haynes et al. 1997; (7) Fisher & Tully 1981; (8) Chamaraux et al. 1990; (9) Peterson 1979; (10) Bicay & 
Giovanelli 1987; (11) Mould et al. 1993; (12) Richter & Huchtmeier 1987; (13) Lewis et al. 1985; (14) Haynes & Giovanelli 1984; 
(15) Williams 1985; (16) Eder et al. 1991; (17) Theureau et al. 2005; (18) Haynes et al. 1999; (19) Bothun et al. 1985; (20) Bicay 
& Giovanelli 1986a; (21) Richter & Huchtmeier 1991; (22) Bicay & Giovanelli 1986b; (23) Williams & Kerr 1981; (24) Hoffman et 
al. 1989; (25) Hoffman et al. 1985: (26) Giovanelli et al. 1997; (27) van Driel et al. 2000; (28) Freudling et al. 1992; (29) Reif et 
al. 1982; (30) Giovanelli & Haynes 1993; (31) Bottinelli et al. 1982; (32) Dean & Davies 1975; (33) Lu et al. 1993; (34) Wegner et 
al. 1993; (35) Nordgren et al. 1998; (36) Huchtmeier 1973: (37) Giovanardi et al. 1983; (38) Bottinelli et al. 1999; (39) Dickel & 
Rood 1978; (40) Magri 1994; (41) Huchtmeier et al. 1995; (42) Giovanardi & Salpeter 1985; (43) Impey et al. 2001; (44) Fosbury 
et al. 1978; (45) Giovanelli & Haynes 1985; (46) Haynes et al. 1988; (47) Staveley-Smith & Davies 1987; (48) Bregman et al. 1992; 
(49) Staveley-Smith & Davies 1988; (50) Bottinelli et al. 1980; (51) Sulentic & Arp 1983; (52) Schneider et al. 1986; (53) Krumm 
& Salpeter 1980; (54) Rosenberg & Schneider 2000; (55) Dell'Antonio et al. 1996; (56) Lewis 1987; (57) Huchtmeier & Richter 
1987; (58) Mirabel & Wilson 1984; (59) Huchtmeier & Richter 1985; (60) Haynes 1981; (61) Huchtmeier 1982; (62) Huchtmeier & 
Seiradakis 1985; (63) Davis & Seaquist 1983; (64) Tifft & Cocke 1988; (65) Balkowski & Chamaraux 1981; (66) Lewis & Davies 
1973; (67) Knapp et al. 1979; (68) Williams & Rood 1987; (69) Bieging & Biermann 1977; (70) Mirabel & Sanders 1988; (71) 
Bosma 1979; (72) Burstein et al. 1987; (73) Haynes & Giovanelli 1991; (74) Schneider et al. 1990; (75) Bottinelli & Gouguenheim 
1977; (76) Balkowski & Chamaraux 1983; (77) Bottinelli & Gouguenheim 1979; (78) Knapp et al. 1985; (79) Freudling 1995; (80) 
Verheijen & Sancisi 2001; (81) Helou et al. 1984; (82) Haynes et al. 1990; (83) Garci'a-Barreto et al. 1994; (84) Kraan-Korteweg 
et al. 1999; (85) Bregman & Roberts 1988; (86) Davies & Lewis 1973; (87) Giovanelli & Haynes 1983; (88) Helou et al. 1981; (89) 
Knapp et al. 1978; (90) van Driel et al. 2000; (91) Hoffman et al. 1987; (92) Krumm & Salpeter 1976; (93) Roth et al. 1994; (94) 
Gavazzi 1987; (95) Rubin et al. 1976; (96) Garcia et al. 1994; (97) Chamaraux et al. 1987; (98) van Driel et al. 2001; (99) Haynes 
& Giovanelli 1991; (100) Staveley-Smith et al. 1992; (101) Roth et al. 1991; (102) Schneider et al. 1992; (103) Chengalur et al. 
1994; (104) Lewis 1983; (105) Barnes et al. 1997; (106) Freudling et al. 1988; (107) Lu et al. 1990; (108) Hewitt et al. 1983; (109) 
Oosterloo & Shostak 1993; (110) Thuan et al. 1999; (111) Irwin & Seaquist 1991; (112) Huchtmeier et al. 1997; (113) Giovanelli & 
Haynes 1985; (114) Matthews et al. 2001; (115) Giovanelli et al. 1981; (116) Thuan & Martin 1981; (117) Schombert et al. 1992; 
(118) Martin et al. 1991; (119) Heckman et al. 1978: (120) Matthews & van Driel 2000; (121) Salzer 1992; (122) Jackson et al. 
1987; (123) Fontanelli 1984; (124) Bothun et al. 1984; (125) Richter & Huchtmeier 1982. 



Hyperleda have been "homogenized" following the pre- 
cepts of McElroy (1995). Generally the measurements 
pertain to a central aperture that is smaller than the effec- 
tive radius of the bulge. Although it is sometimes cus- 
tomary to scale the velocity dispersions to a fixed aper- 
ture (e.g., to /? e ff/8; J0rgensen et al. 1995), this practice 
assumes that all bulges possess a similar velocity dis- 
persion profile, which, according to the observations of 
Pizzella et al. (2005), appears not to be the case. We 
thus use the central values of uq, which for our sample 
on average have uncertainties of 13%, with a standard 



deviation of 1 1%. 

Nearly all of the rotation velocities in Hyperleda are de- 
rived from spatially integrated H I line widths, homoge- 
nized in the manner described in Paturel et al. (2003b). 
Using extensive sets of galaxies that contain both inte- 
grated H I profiles and spatially resolved optical rotation 
curves, these authors determined the optimal transfor- 
mation between the H I line width measured at different 
levels of the line peak (e.g., 20% or 50% of the max- 
imum) and the maximum velocity of rotation. In this 
way, the archived H I measurements, which can be quite 
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voluminous for any given galaxy, can be standardized 
and reduced to a single value for the rotation velocity. 
Although this procedure is clearly very useful, following 
the spirit in which the stellar velocity dispersions them- 
selves were treated, unfortunately the rotation velocities 
listed in Hyperleda cannot be used at face value because 
of the possibility of source confusion. Single-dish H I 
measurements subtend over a significant beam area, typ- 
ically from 3. '5 for the Arecibo telescope to as much as 
22' for the Nancay telescope, within which there is a 
nonnegligible probability of contamination from inter- 
lopers or neighboring galaxies. Hyperleda does not take 
this crucial effect into consideration prior to homogeniz- 
ing the H I data. 

Accordingly, we have taken the following steps to cull 
the data. Beginning with an initial sample of ~ 1500 
galaxies having both velocity dispersions and H I mea- 
surements, for each galaxy we systematically inspected 
the list of H I line widths given in Hyperleda. In this pro- 
cess, we give preference to more modern observations if 
available, to data taken with the highest spectral resolu- 
tion, to line width measurements that pertain to either the 
20% or 50% of the maximum of the line profile (prefer- 
ably both), and, whenever possible, we try to minimize 
the heterogeneity of the sample by favoring larger, more 
systematic surveys. For each galaxy, we identify what 
we deem to be the most robust line width measurement 
taken with the smallest available beam, where robust- 
ness is judged by whether the line width has reached a 
stable, asymptotic value when data from multiple tele- 
scopes are listed. We then carefully examined digital 
optical images, from the SDSS if available or else from 
the scanned images of the Palomar Digital Sky Survey, 
in combination with redshift information listed in the 
NASA/IPAC Extragalactic Database (NEDfl to elimi- 
nate potential sources of confusion within a radius equal 
to the full-width at half power of the beam of the cho- 
sen telescope (see Table 2). For observations taken with 
the Arecibo telescope, the search radius was increased 
to 7f5 to account for the extended sidelobes of the beam; 
the intensity of the first sidelobes of the beam drops to 
~10% of the peak at a distance of 5. '5 from the beam 
center, and by 7'-8' it becomes negligible (Heiles et al. 
2000). A number of galaxies contain faint, low-surface 
brightness companions that formally lie within the H I 
beam, but generally such companions can be ruled out 
as sources of confusion because of their low luminosities 
(they would otherwise grossly violate the Tully-Fisher 
relation given the large line widths). Galaxies that are 
likely to be confused are flagged and omitted from the 
sample. We retained a few galaxies that, although po- 
tentially confused within the single-dish beam, have line 
widths consistent with the velocity amplitude measured 
in interferometer maps, which suggest that confusion is 
not a problem. To avoid complications in the interpre- 
tation of their H I kinematics, we also removed galax- 
ies known to be merger remnants or that, from our vi- 
sual inspection, otherwise possess disturbed morpholo- 
gies, tidal tails, shells, or polar rings. We do not consider 
dwarf irregular galaxies (T > 9.0) because they are un- 
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likely to have trustworthy central stellar velocity disper- 
sion measurements or unambiguous bulges. The above 
selection reduced the Hyperleda sample to 293 galaxies. 

• Sample 2 — Many galaxies in Hyperleda have usable 
H I line widths but no stellar velocity dispersion in the 
database. A significant number of them overlap with 
the Sloan Digital Sky Survey (SDSS; York et al. 2000), 
and stellar velocity dispersions for 435 of these have 
been measured and made publically available by D. J. 
Schlege0. In this selection, we only retain galaxies that 
have velocity dispersions larger than the spectral instru- 
mental resolution of SDSS (~ 70 km s" 1 ) and that satisfy 
our morphological type cut (T < 9.0). We visually ex- 
amined the SDSS spectrum of every object to confirm 
that the reported velocity dispersion is sensible; a num- 
ber of spurious values, almost all resulting from erro- 
neous fitting of noisy spectra of faint, dwarf galaxies, 
were rejected during this process. D. J. Schlegel et al. 
(2007, in preparation; see also Heckman et al. 2004) 
measure velocity dispersions using a direct-fitting algo- 
rithm and template stars from a library of echelle stel- 
lar spectra. The Appendix presents an external compari- 
son of the velocity dispersions given in SDSS with those 
published in Hyperleda, for objects in common between 
samples 1 and 2. The SDSS values are on average 12% 
larger than those in Hyperleda. 

• Sample 3 — The issue of whether ao or v m more funda- 
mentally tracks Mbh can be most effectively addressed 
by using extremely late-type spirals that essentially lack 
a bulge altogether (Boker et al. 2003) and rarely con- 
tain a central black hole (see discussion in §4.1). While 
rotational velocities can be measured readily in these 
gas-rich systems, obtaining accurate central stellar ve- 
locity dispersions for them is extremely challenging be- 
cause of the faintness or absence of a central spheroidal 
component. We can place a useful limit on the "bulge" 
velocity dispersion of such systems by measuring the 
dispersion of its central nuclear star cluster, which is 
commonly found (Boker et al. 2002). These are chal- 
lenging observations, because the clusters are faint and 
the dispersions are small, requiring echelle resolution on 
large telescopes. We have been able to locate suitable 
Co measurements for a total of 10 late-type, bulgeless 
spirals: NGC 598 (M33; Kormendy & McClure 1993), 
NGC 4395 (Filippenko & Ho 2003), and eight galax- 
ies from the Boker et al. (2002) survey (Walcher et al. 
2005). 

• Sample 4 — Courteau et al. (2007) give the latest com- 
pilation of galaxies having both Co and v,„ determined 
largely from extended optical rotation curves [59 of their 
objects come from the study of Prugniel et al. (2001), 
whose measurements of v m were derived from H I line 
widths]. Of the 164 galaxies in Courteau's compilation, 
54 are not included in our samples 1-3; these form the 
last sample in our study. We distinguish three subsets: 
(a) 21 galaxies with ctq taken from Courteau et al. and 

v m derived from H I line widths in Hyperleda; (b) 26 




galaxies with v m taken from Courteau et al. and <jq 
taken from Hyperleda; and (c) 7 galaxies with <jq and 
v m both taken from Courteau et al. In the interest of 
maximizing the homogeneity of our combined sample, 
we have decided to use H I line widths whenever avail- 
able even if Courteau et al. give v m from resolved rota- 
tion curves. Similarly, we give preference to the ao val- 
ues in Hyperleda because these represent weighted av- 
erages of all available literature measurements, and they 
have been scaled to a common system adopted for the 
rest of our study. The velocity dispersions that overlap 
between Courteau and Hyperleda show excellent agree- 
ment, with Courteau's values on average 3% higher than 
Hyperleda' s. 

The combination of the above four samples produces a final 
compilation of 792 galaxies, summarized in Figure 1 and Ta- 
bles 1 and 2. Most of the galaxies are nearby (median distance 
23 Mpc) and luminous (median M B \ = -20.3 mag), spanning 
the entire range of Hubble types, from giant ellipticals to late- 
type spirals (morphological type index T = -5 to 8). 

The published H I line widths have been corrected for in- 
strumental resolution but not for (1 + z) cosmological stretch- 
ing or for broadening by interstellar turbulence. A number 
of authors have discussed the impact of turbulence on the line 
width measurements (e.g., Tully & Fouque 1985; Fouque et al. 
1990). Given the demographic make up of most of our sample 
(Fig. lb), we adopt the simple linear subtraction recommended 
by Bottinelli et al. (1983) for giant galaxies, 

W = W obs -W turb , (1) 

where W b s is the observed line width (at either 20% or 50% 
of the maximum), and the turbulent velocity is taken to be 
W turb = 22 km s" 1 for W 2 q and W turb = 5 km s" 1 for W 50 (Ver- 
heijen & Sancisi 2001). 

As the H I line widths must be deprojected along the line- 
of-sight, we need to pay careful attention to the adopted in- 
clination angle. Although Hyperleda conveniently lists inclina- 



tion angles for our sample, they differ systematically and some- 
times dramatically from the values given in other catalogs (see 
Appendix). We therefore resorted to recompiling our own ax- 
ial ratios and recalculating the inclination angles for the entire 
sample (Table 1). We derive the inclination angle i using Hub- 
ble's (1926) formula, 



l ~1o 

which makes use of the apparent flattening of the galaxy (q = 
b/a) as measured from the ratio of its semi-minor and semi- 
major isophotal diameters at a surface brightness of /z^ = 25 
mag arcsec" 2 , assuming that the intrinsic thickness of the disk 
(go) depends on morphological type as given in Paturel et al. 
(1997). To mitigate against large uncertainties inherent in 
determining the inclination angle of nearly face-on systems, 
throughout this paper we remove galaxies that have inclination 
angles of i < 30°; these account for <~ 10% of the original sam- 
ple. An important caveat, however, is that for most of the galax- 
ies in our study we have absolutely no information on the spatial 
distribution of the neutral hydrogen. Given that some H I disks 
can be misaligned with respect to the stellar distribution (e.g., 
Lewis 1987), our optically derived inclination correction may 
not be valid in some instances. Another source of uncertainty 
comes from strong warps and other types of nonaxisymmetric 
distortions in the H I disk, which are known to be present in 
many disk galaxies (e.g., Baldwin et al. 1980; Bosma 1981; 
Lewis 1987; Richter & Sancisi 1994; Haynes et al. 1998). In- 
deed, in §4.3 we present evidence based on our sample that a 
sizable fraction of nearby galaxies may have strongly disturbed 
or dynamically unrelaxed H I distributions. 

Lastly, we need an empirical prescription to convert the cor- 
rected H I line width into the maximum rotation velocity. This 
issue has been investigated in a number of studies (e.g., Math- 
ewson et al. 1992; Courteau 1997), the latest and most thor- 
ough being that of Paturel et al. (2003b), whose calibrations we 
adopt. In detail, the calibrations depend on the resolution of the 



6 



Bulge and Halo Kinematics 




FIG. 2. — The relationship between central stellar velocity dispersion and the maximum rotation velocity for galaxies of the following Hubble types: (a) all, (b) 
ellipticals (-5 <T< -2.5), (c) SO (-2.5 < T < -0.5), (d) Sa-Sb (-0.5 < T < 2.5), (e) Sbc-Sc (2.5 <T< 4.5), (/) later than Sc (T > 4.5), (g) all disks, and (h) all 
spirals. Sources that are low-velocity and high-velocity outliers in the Tully-Fisher relation (Fig. 3) are plotted as open triangles and squares, respectively. The solid 
points belong to the "kinematically normal" sample. In panels (a) and (f)-(h), the extreme late-type, bulgeless spirals are plotted as stars. For the sake of clarity, 
error bars are not shown in panels (a), (g), and (h). The dashed lines denote, from bottom to top, v m = uq, v2 uq, and \/4 <ro. The best fit to the entire sample of 
kinematically normal galaxies is plotted as a solid line in panel (a). 



observations, but for simplicity we just use those obtained for 
the highest resolution (8 km s" 1 ): 

log 2w ffl sin; = (1.071 ± 0.009) logW 50 - (0.210 ± 0.023) 

log2u ffl sin/ = (1.187±0.002)logW 20 -(0.543±0.005). (3) 

The uncertainty on v m is taken to be the value formally given 
in Hyperleda. This is a conservative estimate because in many 
cases the Hyperleda value incorporates a larger spread of indi- 
vidual measurements than we actually used. For our sample, 
the values of v m have an average uncertainty of 6.3%, with a 
standard deviation of 5.4%. 

3. THE v m -<ro RELATION 

Figure 2 summarizes our principal findings. Taken collec- 
tively (Fig. 2a), the distribution of v m versus Co shows, at best, 
a loose correlation. Although the best-fit relations of Baes et 
al. (2003) or Pizzella et al. (2005) roughly bisect the cloud of 
points, the scatter is enormous, far greater than can be attributed 
to observational errors or potential sources of systematic uncer- 
tainty (e.g., inclination corrections for v m or aperture correc- 
tions for cto). Separating the sample by Hubble type (Fig. 2b- 
2f) reveals three main culprits for the large scatter: (1) a system- 
atic shift of zeropoint as a function of Hubble type, most clearly 
seen in the locus of the ridgeline defining the upper envelope of 
the distribution of points; (2) the existence of subset of low-oo 
(cto <, 50 km s" 1 ), extreme late-type galaxies (r w 6-9, Hub- 
ble types Scd-Sdm) that have almost constant v m (~ 100 km 



s ); and (3) a cloud of points, present in all Hubble type bins, 
but especially prominent for earlier-type systems, with very low 
values of u m /(To (plotted as open triangles). 

To assess how much of the scatter in the v m -ao diagram is 
intrinsic, we make use of the Tully-Fisher relation to constrain 
what value of v m any particular galaxy ought to have given its 
luminosity. Since the Hubble type mix of a sample affects the 
slope, normalization, and scatter of the Tully-Fisher relation 
(e.g., Roberts 1978; Rubin et al. 1985; Bell & de Jong 2001; 
De Rijcke et al. 2007; Pizagno et al. 2007), and the variations 
are minimized in the near-infrared (e.g., Verheijen 2001), we 
assembled /T,-band (2.16 /im) magnitudes for nearly the entire 
sample using the Extended Source Catalog of the Two-Micron 
All-Sky Survey (2MASS; Skrutskie et al. 2006). The photome- 
try pertains to the "total" magnitudes (Table 1), and for simplic- 
ity we do not correct for Galactic or internal extinction, which 
should be quite small in the K s band. Figure 3 shows that a K s - 
band Tully-Fisher relation exists for all Hubble types, including 
E and SO, and, importantly in the present context, that there is 
little obvious variation in the slope, normalization, or scatter 
of the relation across the wide range of Hubble types included 
in our sample. Overplotted on the figure is the fit derived by 
Verheijen (2001) for the K 1 band, which is quite similar to the 
2MASS K s band (Bessel 2005). While the observed scatter of 
our Tully-Fisher relation is larger than that found by Verheijen 
(2001), a result that can be anticipated considering the larger 
distance errors and greater heterogeneity of our sample, most 
of our objects (625/792 or 70%; solid points) fall comfortably 
within the boundaries that enclose twice the rms scatter 



Ho 



7 




log (v m / km s ! ) log {v m / km s l ) log (v m / km s : ) log {v m / km s" 1 ) 



FIG. 3.— Tully-Fisher relation in the K s band for galaxies of the following Hubble types: (a) all, (b) ellipticals (-5 < T < -2.5), (c) SO (-2.5 < T < -0.5), (d) 
Sa-Sb (-0.5 < T < 2.5), (e) Sbc-Sc (2.5 <T< 4.5), (/) later than Sc (T > 4.5), (g) all disks, and (h) all spirals. The solid line represents the £'-band Tully-Fisher 
relation from Verheijen (2001); the dashed lines mark the region that has twice the rms scatter. Sources that are low- velocity and high- velocity outliers are plotted 
as open triangles and squares, respectively; the solid points are considered to be the "kinematically normal" members of the sample. In panels (a) and (f)-(h), the 
extreme late-type, bulgeless spirals are plotted as stars. For the sake of clarity, error bars are not shown in panels (a), (g), and (h). 



of Verheijen's "H I" sample (0.59 mag); hereafter we will re- 
fer to these as the "kinematically normal" objects. Neverthe- 
less, Figure 3 shows that there are a significant number of out- 
liers, which we loosely and somewhat arbitrarily define to be 
those that lie outside of the 2-rms band. Closer inspection re- 
veals that there is an excess of low-velocity objects (132/792 or 
17%; open triangles) compared to high-velocity ones (35/792 
or 4%; open squares). This is most obvious among the spirals 
(T > -0.5) as a low-w m "plume," but to a lesser extent it can be 
seen also among the E and SO systems. 

Within this backdrop, the distribution of points in the v m -ao 
diagram can now be more easily interpreted. The "kinemati- 
cally normal" objects on the Tully-Fisher relation obey a loose 
correlation between v m and ao, roughly occupying the region 
from v m = ao to v m = Viao, depending on Hubble type. The 
correlation has significant scatter, and it is nonlinear. An or- 
dinary least-squares bisector fit for the entire sample (Fig. 2a) 
yields 

logu ffl = (0.82 ± 0.027) logcr + (0.57 ±0.058), (4) 

very similar to the fit reported by Ferrarese (2002). There is no 
significant variation of the slope with Hubble type. Limiting 
the fit to the 550 kinematically normal spiral galaxies, the fit is 
nearly identical: 

log v m = (0.80 ± 0.029) logcr + (0.62 ± 0.062). (5) 

The bulgeless late-type galaxies, which are largely kinemati- 
cally normal in the Tully-Fisher relation, continue to depart no- 
tably from the rest of the sample in the v m -<Jo diagram. It is 



of interest to note that the few late-type spirals (NGC 6140, 
NGC 6689, PGC 28990) that do host small bulges (as op- 
posed to nuclear star clusters), along with some dwarf ellipti- 
cals and SOs (NGC 3870, NGC 7077, PGC 5441, PGC 71938), 
still follow the low-velocity extrapolation of the v m -ao relation. 
Within the large scatter, we find no compelling evidence for a 
break in the v m -ao relation at low velocities. Objects with low 
v m lo~o (<;1) comprise almost exclusively the low-velocity out- 
liers in the Tully-Fisher relation. The nature of the \cw-v m lao 
objects is discussed further in §4.3. The high-velocity outliers, 
on the other hand, occupy the upper envelope of the v m -ao rela- 
tion, especially for ao ^> 80 km s _1 . Had these objects not been 
excluded, they would bias the v m -ao relation to a steeper slope. 

Figure 4 highlights the trends with Hubble type more explic- 
itly, showing the distribution of vjao for six bins of morpho- 
logical types. Although the binning of the morphological types 
is somewhat arbitrary, nevertheless the pattern is clear: as the 
Hubble type becomes later, the median of the distribution of 
v m /ao systematically shifts to larger values, from ~1.2 for el- 
lipticals, to <~1.4 for early-type spirals, to ~1.8 for late-type 
spirals, and, very dramatically, to ^4.6 for extreme late-type, 
bulgeless spirals. This trend was already largely noticed by ear- 
lier studies (Whitmore et al. 1979; Whitmore & Kirshner 1981; 
Franx 1993; Whittle 1992a; Zasov et al. 2005). Since we do 
not have reliable bulge-to-disk photometric decompositions for 
the majority of our sample, we constructed a surrogate measure 
of bulge luminosity for the disk galaxies by using the observed 
values of ao in combination with the Faber-Jackson (1976) re- 
lation. Using the 1072 elliptical (T < -2.5) galaxies with mea- 
surements of ao listed in Hyperleda, we find the following 
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FIG. 4. — The distribution of the ratio of disk rotation velocity to bulge ve- 
locity dispersion for the "kinematically normal" galaxies (see §3) of Hubble 
type E, SO, Sa-Sb, Sbc-Sc, Scd-Sdm, and extreme late-type, bulgeless spirals. 
The median of each distribution is given, followed by the number of galaxies 
in the group; the median is also marked by an arrow. 

ordinary least-squares bisector fijj Mb t = -6.801ogcro-4.89. 
This method of estimating the "bulge" luminosity certainly has 
limitations. The Faber-Jackson relation for bulges is, as yet, not 
well-determined: it appears to differ systematically with Hub- 
ble type (Whitmore & Kirshner 1981; Kormendy & 
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FIG. 5. — The variation of the ratio of disk rotation velocity to bulge velocity 
dispersion as a function of the bulge-to-total luminosity ratio. Here the "bulge" 
luminosity is calculated from the Faber-Jackson relation, as determined from 
a sample of 1072 elliptical galaxies in Hyperleda. The extreme late-type, bul- 
geless spirals are plotted as stars. The solid points connected by the dashed 
line show the mean and standard deviation of v m /ao binned by B/T for the 
"kinematically normal" sources (see §3). 



Illingworth 1983; but see Whittle 1992b) as well as between 
classical bulges and "pseudo-bulges" (Kormendy & Kenni- 
cutt 2004), and it shows significant intrinsic scatter within 
each Hubble type. Nevertheless, it is debatable whether these 
sources of systematic uncertainty are any more serious than 
the vagaries of photometric bulge-to-disk decomposition, es- 
pecially for such a large and diverse sample of galaxies. With 
these caveats in mind, Figure 5 shows that v m /(To clearly varies 
systematically with bulge-to-total (B/T) luminosity raticQ. A 
given value of B/T can host a large range of v m /<ro, but on 
average v m /ao increases as the prominence of the bulge com- 
ponent decreases, qualitatively consistent with the trends noted 
in Figures 2 and 4. 

Another crude, but more direct, method to estimate the de- 
gree of bulge dominance is to calculate the "concentration in- 
dex" of the light profile, which is roughly related to the bulge- 
to-disk ratio or Hubble type (e.g., Doi et al. 1993). The 
correlation between concentration index and Hubble type in 
the SDSS database has been investigated by Shimasaku et al. 
(2001) and Strateva et al. (2001). This is the approach taken 
by Courteau et al. (2007), who defined a concentration index 
C28 = 5 log(r8o / 7"2o), where r2o and r%o are the radii that enclose 
20% and 80% of the total light, respectively. They computed 
C28 for 81 out of the 164 galaxies in their sample with available 
/-band images in the SDSS. We follow Courteau et al. in using 
the concentration index as a surrogate indicator for the den- 
sity distribution of the galaxy light profile, but instead of C28, 
we simply adopt the cataloged /-band Petrosian radii enclosing 
50% (rpso) and 90% (rp9o) of the light to form an equivalent 
concentration index C59 = 51og(rp9o/rpso)- Not surprisingly, 
C28 and C59 are strongly correlated (Fig. 6a), albeit with signifi- 
cant scatter, which, too, is not unexpected. Figure 6b illustrates 
that substituting C59 for C28 faithfully recovers the correlation 
between v m /ao and concentration index presented in Courteau 
et al. (2007; see their Fig. 2). If anything, using C59 instead 
of C28 seems to produce a diagram with even somewhat less 
scattefl 

The majority of our sample (76%) have photometric data 
available in the Fifth Data Release of SDSS (Adelman- 
McCarthy et al. 2007), and so we can calculate the C59 parame- 
ter for these objects (Table 1). As anticipated from the Courteau 
et al. study, Figure 7 illustrates that our sample indeed also be- 
haves very similarly: a systematic trend exists between v m /ao 
and C59, although the scatter is substantial. 

The density distribution, as traced by the morphological type, 
bulge-to-disk ratio, or concentration index, all of which are 
loosely mutually correlated, appears to be the dominant factor 
that determines a galaxy's v m /ao ratio. We have examined the 
total galaxy luminosity as a possible additional parameter, but 
it has little or no effect. This is not surprising, because with the 
exception of very late-type systems, a galaxy's total luminosity 
shows only subtle variation with Hubble type (e.g., Roberts & 
Haynes 1994). A galaxy's integrated broad-band colors (U-B, 
B-V, B-I, B- K s ), to the extent that they are available from 
Hyperleda, do show a moderate correlation with v,„/ao, but this 
is likely just a reflection of the dependence of galaxy color on 



4 This fit is roughly consistent with the SDSS results of Bernardi et al. (2003). For example, for a velocity dispersion of <ro = 200 km s 1 , their equation (11) pre- 
dicts M g = —19.33 mag, which, for g-B = -0.53 mag expected for elliptical galaxies (Fukugita et al. 1995), translates to Mb = -18.80 mag. Our fit gives Mb = —18.49 
mag. 

5 A small fraction of the sample, not shown in Figure 5, has B/T values that formally exceed unity. This should not be viewed as too alarming, in view of the 
indirect method by which we have estimated the "bulge" luminosity. 

6 Figure 6b contains fewer data points than Figure 2 of Courteau et al. (2007) because the latter treated repeat observations of the same objects as independent 
galaxies (S. Courteau 2006, private communications). 
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FIG. 6. — (a) The relation between i-band concentration indices as derived 
in the Petrosian system of SDSS, C59, and as given by Courteau et al. (2007), 
C28- The two concentration indices are strongly correlated, (b) The variation 
of v m /o-Q with C59 for the sample of Courteau et al., plotted in the same fashion 
as Fig. 2 of Courteau et al., demonstrating that C59 can be substituted for Ci%. 

morphological type or concentration. 

One caveat should be kept in mind. It is possible that the 
central velocity dispersions of the later-type galaxies are under- 
estimated due to contamination by disk light. If so, this could 
mimic the variation of v m /<7o with galaxy type or bulge-to-disk 
ratio. This possibility needs to be checked with detailed spa- 
tially resolved stellar kinematical data. 

4. DISCUSSION 

4. 1 . Implications for Black Hole Demographics 

This study was initiated with the intention of better quantify- 
ing the relation between disk rotation velocity and bulge stellar 
velocity dispersion, with the hope that such a correlation might 
be used as a new empirical tool to study the demographics of 
central black holes in galaxies. Earlier studies (Ferrarese 2002; 
Baes et al. 2003), based on relatively small samples with ex- 
tended optical rotation curves, found that galaxies over a wide 
range of Hubble types, from mid-type spirals to giant ellipti- 
cals, obey a strong correlation between v m and 00 with almost 
no intrinsic scatter, at least for v m ;> 80 km s" 1 . Subsequent 
work has emphasized that the v m -ao relation may depend on 
surface brightness, with low-surface brightness galaxies lying 
systematically offset (larger v m for a given 00) from a still tight 
sequence defined by high-surface brightness galaxies (Pizzella 
et al. 2005; Buyle et al. 2006). If high-surface brightness galax- 
ies truly do obey a tight v„-(jq relation, one would have to 
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FIG. 7. — The variation of u m /o"o with the concentration index C59. The 
extreme late-type, bulgeless spirals are plotted as stars. The solid points con- 
nected by the dashed line show the mean and standard deviation of v m /ao 
binned by C59, for the "kinematically normal" sources (see §3). 

reevaluate which dynamical variable (cto or v m ) or, equivalently, 
which structural component of the galaxy (bulge or halo) is 
more fundamentally related to black hole mass. 

Nevertheless, there is reason to suspect that the v m -ao rela- 
tion may not be as simple as had been claimed. Indeed, from 
the initial work of Whitmore and collaborators (Whitmore et al. 
1979; Whitmore & Kirshner 1981), as well as from subsequent 
updates thereof (Fall 1987; Whittle 1992a; Franx 1993; Nelson 
& Whittle 1996), one could already have concluded that the 
ratio v m /ao is not a constant, but rather that it varies systemati- 
cally as a function of bulge-to-disk ratio, in the sense that early- 
type galaxies have characteristically smaller vJctq ratios than 
late-type galaxies. It was also apparent that there is significant 
scatter in v m /ao for any given bulge-to-disk ratio. Baes et al. 
(2004) speculated that the culprit for the increased scatter could 
lie in Whitmore's use of integrated H I line widths, which may 
not accurately trace the flat part of the rotation curve. This ex- 
planation, however, seems implausible in light of the long his- 
tory of Tully-Fisher studies that have made use of integrated H I 
line widths, not to mention of the many comparisons made with 
optically derived rotation curves (e.g., Rubin et al. 1978; Thon- 
nard 1983; Mathewson et al. 1992; Courteau 1997; Paturel et 
al. 2003b). Our own comparison shown in Figure 9, though 
limited, reinforces previous conclusions that integrated H I line 
widths serve as an effective proxy for rotation velocities mea- 
sured from extended rotation curves, at least for spiral galaxies. 
Moreover, even if one were to admit that H I-based rotation ve- 
locities were less accurate, one would still have to explain why 
v m /ao changes systematically with bulge-to-disk ratio. Zasov 
et al. (2005) provided further evidence that the v m -ao rela- 
tion contains significant intrinsic scatter. Using a collection of 
41 galaxies with black hole masses, central stellar velocity dis- 
persions, and rotational velocities from optical rotation curves, 
these authors concluded that although Mbh does correlate with 
v m , it is better correlated with <to. They also find that, at a 
given v m , SO-Sab galaxies tend to have larger M B h than later- 
type galaxies. From this, one can infer that earlier-type galaxies 
possess a smaller v m /o-Q ratio than later-type galaxies. The lat- 
est demonstration that v m /ao varies considerably and system- 
atically with galaxy type comes from Courteau et al. (2007), 
whose v m compilation — now increased to 164 galaxies — are 
also based largely on spatially resolved optical rotation curves. 
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Their study shows that the v m -ao relation exhibits marked scat- 
ter and that v m lao decreases steadily, albeit with large fluctua- 
tions, with increasing galaxy light concentration. 

The present study evaluates the v m -cro relation with the 
largest sample to date, totaling 792 galaxies. This is made pos- 
sible not only by access to a large number of new cto measure- 
ments from SDSS, but also by taking advantage of integrated 
H I line widths to estimate rotation velocities, a shortcut whose 
efficacy has been amply proved. As anticipated, we find that 
the v m -cro relation is far from tight. The zeropoint shifts sys- 
tematically as a function of Hubble type, and within each Hub- 
ble type bin the scatter is at least a factor of 2-3. From the 
point of view of using v m to predict Mbh, this is not promising, 
not if other properties of the host galaxy such as bulge lumi- 
nosity (Kormendy & Gebhardt 2001; Marconi & Hunt 2003), 
and especially ao (Tremaine et al. 2002), are available. On the 
other hand, there are instances when the stellar component of 
the host galaxy can be extraordinarily challenging, if not impos- 
sible, to detect, either kinematically or photometrically. This 
is often the case in luminous active galactic nuclei or quasars. 
The host galaxy's neutral hydrogen gas, by contrast, being im- 
mune to the bright glare of the active nucleus, may be more 
readily detectable with deep radio observations. Under these 
circumstances, an integrated line width, which can be extracted 
from an H I spectrum of even moderate signal-to-noise ratio, 
may be the only constraint available on the host galaxy. The 
recent analysis by Ho et al. (2007b) of a new H I survey of 
active galaxies (Ho et al. 2007a) demonstrates the rich body 
of information on the host galaxy that can be ascertained in 
this unique fashion. The existence of low-w„,/<7o outliers (§4.3) 
certainly complicates matters, but this difficulty is not insur- 
mountable because these objects can be recognized as outliers 
in the Tully-Fisher relation, if an independent estimate of the to- 
tal host galaxy luminosity can be estimated. Alternatively, Ho 
et al. (2007b) suggest that this kinematically anomalous pop- 
ulation may have a greater tendency to exhibit single-peaked 
and/or highly asymmetric line profiles. While this trend still 
needs to be verified with better data, it may offer a very effec- 
tive, practical strategy to weed out unwanted contaminants. 

It is important to stress that the large scatter of the v m -ao 
relation is not a consequence of mixing low- and high-surface 
brightness galaxies. Although we cannot rigorously identify 
low-surface brightness galaxies in our sample with the material 
at hand, given the very bright luminosities of most of our sam- 
ple it is highly improbable that a large fraction of them are low- 
surface brightness galaxies (see, e.g., Sprayberry et al. 1997). 

The absence of a single, universal u,„-ero correlation, in con- 
cert with the general consensus that the Mbh-oo relation is 
tight, casts doubt on Ferrarese's (2002) hypothesis that the halo 
mass is more fundamentally connected to the black hole mass 
than is the bulge mass. Judging from both direct dynamical 
searches of black holes (e.g., Magorrian et al. 1998; Kormendy 
2004, and references therein) and the high detection rate of ac- 
tive nuclei in early-type galaxies (Ho et al. 1997; Ho 2004), the 
occupation fraction of black holes in bulges is very high, per- 
haps approaching unity. By contrast, the inverse may be true 
in late-type and dwarf galaxies. None has so far yielded a di- 
rect dynamical detection of a central black hole (Gebhardt et 
al. 2001; Merritt et al. 2001; Valluri et al. 2005), and nu- 



clear activity is exceedingly rare (Ho et al. 1997; Ulvestad & 
Ho 2002), though not unknowrQ (Filippenko & Ho 2003; Barth 
et al. 2004, 2005; Greene & Ho 2004, 2007a, 2007b). Thus, 
with a few exceptions, central black holes, be they active or 
inactive, are associated almost exclusively with bulges, not ha- 
los. The very late-type spirals included in this study (sample 3) 
possess substantial rotation velocities (65 < v m < 140 km s" 1 ) 
and correspondingly non-negligible halo masses but effectively 
no bulges (Boker et al. 2003): their central stellar velocity dis- 
persion, 13 < (To < 34 km s" 1 , arises not from a bulge but a 
compact, nuclear star cluster (Kormendy & McClure 1993; Fil- 
ippenko & Ho 2003; Walcher et al. 2005). If M BH were more 
fundamentally linked with v m than do, we would expect these 
galaxies to possess black holes with Mbh ~ 2 x 10 5 to 4 x 10 6 
Mq [using the v m -ao relation of Baes et al. (2003) and the 
Mbh-oo relation of Tremaine et al. (2002)], and yet none but 
NGC 4395 is known for sure to contain a central black hole (Fil- 
ippenko & Ho 2003). M33 (NGC 598) has the most stringent 
limit: M BH < 1500 M Q (Gebhardt et al. 2001). 



4.2. Implications for Galaxy Formation 

Discouraging as it may be as a black hole mass predictor, the 
v m -(Jo relation can be regarded as a useful, additional constraint 
on galaxy formation models. What factors determine the v m /ao 
ratio in any given galaxy? Our study, in agreement with that 
of Courteau et al. (2007), shows that v m /ao depends, at least 
in part, on the mass distribution within a galaxy. This is evi- 
dent in the correlation between v m /ao and three mutually related 
quantities that reflect the density profile of a galaxy, namely its 
Hubble type, bulge-to-disk ratio, and concentration index. The 
trends are systematic and statistically significant, but none of 
the correlations can be regarded as particularly tight, suggest- 
ing that v m lo-Q probably depends on more than just the density 
profile alone. This can be anticipated from consideration of the 
collisionless Boltzmann equation for a stellar system in gravi- 
tational equilibrium, in which the rotation velocity and velocity 
dispersion at radius r are related by 



v 2 (r) = aj(r) 



d\np(r) d\na^(r) 



d\nr 



d\nr 



■2 1 



A \ 1 



(6) 



Here a r and ag are the velocity dispersions along the radial and 
tangential directions, and p is the stellar density profile. Fol- 
lowing Courteau et al. (2007), we can recast this equation as 



i-'O 



dlnp(r) d\Ti(j}.(r) 



d\nr 



d\nr 



■2 1 



1/2 



(7) 



where vq denotes the rotation velocity in the inner part of the 
galaxy. The dependence of v m /ao on the density profile enters 
directly through the first term in the brackets as well as through 
v m /vo, which depends on the shape of the rotation curve and 
hence on the detailed mass distribution of the galaxy (Rubin et 
al. 1978, 1985). 

The challenge for galaxy formation models is not to ex- 
plain why v m /ao varies (as it is obvious from equation 7 that it 
should), but rather why galaxies in the local Universe populate 



the v m -ao plane in the manner observed. In this sense, we can 

'Previous authors (Ferrarese 2002; Baes et al. 2003; Buyle et al. 2006) have used the apparent break down of the v,„-ao relation at low velocities to infer that 
low-mass galaxies do not host central black holes. This is incorrect. Bona fide examples of active galactic nuclei with low-mass black holes in low-mass (even dwarf) 
galaxies are known. Although such objects are apparently rare, current search techniques introduce strong selection biases and preclude any definitive conclusions 
regarding their true space densities (Greene & Ho 2007a). 



Ho 



11 



regard the distribution of galaxies on the v m -ao diagram as an- 
other local empirical boundary condition — much like other em- 
pirical galaxy scaling relations — that can be used to both guide 
and constrain theoretical models. Courteau et al. (2007), for ex- 
ample, attempted to reproduce, with limited success, the actual 
variation of vjo^ with concentration index using the equilib- 
rium disk-bulge-halo models of Widrow & Dubinski (2005). 

4.3. Origin of the Low-v m /ao Systems 

One of the most notable features in the v m -ao diagram is 
the existence of a significant number of outliers with very low 
values of v m /ao (open triangles in Fig. 2). As an operational 
definition, we have designated this population as those that de- 
viate from the /f,-band Tully-Fisher relation toward the low- 
v m side by more than twice the rms scatter. This definition, 
though somewhat ad hoc, effectively isolates most of the mem- 
bers characterized by vJuq <; 1 . Although the low-u m /cro ob- 
jects comprise only 17% (132/792) of the total sample, they 
are nonetheless very intriguing because the kinematic measure- 
ments are robust. Some of these systems have rotation veloc- 
ities as low as ^30 km s" 1 , and yet most are luminous, with 
M& ~ -19 to -21 mag. Moreover, these galaxies have central 
stellar velocity dispersions as large as cto ~ 100-250 km s" 1 : 
it is impossible for a cold disk rotating this slowly to survive in 
such a dynamically hot system. The low-t> m /oo outliers can be 
found in essentially all Hubble type bins included in our study, 
but they appear to be somewhat more prevalent among early- 
type (E and SO) galaxies: they make up 33% of the ellipticals, 
40% of the SOs, and ~ 15% to 25% of the spirals, depending 
on Hubble type. 

Apart from their low vJvq ratios, the outliers do not stand 
out in any other obvious manner, at least with the data we read- 
ily have at our disposal (using data given in Tables 1 and 2 
and as listed in Hyperleda). When normalized with respect to 
other galaxies of the same morphological type, they appear to 
be normal in terms of absolute magnitude, mean surface bright- 
ness (measured either at the effective radius or within /ig = 25 
mag arcsec" 2 ), integrated broad-band colors (U-B, B-V, B-I, 
B-K s ), and H I content (either the total H I mass, M H i, or the H I 
mass normalized to the B-band luminosity, Mu JLb). Among 
the disk galaxies in the sample, 56% are classified as barred, to 
be compared with 42% for the kinematically normal sample; it 
is unclear if the marginal excess of barred galaxies in the for- 
mer sample is statistically significant. Other than the indirect 
effect of the morphology-density relation (e.g., Dressier 1980), 
no obvious trends with environment exist, although except for 
noting cluster membership we did not attempt to quantify the 
environment in any rigorous manner. 

How might a luminous galaxy with a sizable hot spheroidal 
component have such kinematically cold neutral hydrogen? 
The first and most obvious possibility is that we have severely 
misjudged the inclination correction to the H I line width. Re- 
call that the inclinations were derived from the axial ratios of 
the stellar isophotes, which implicitly assumes that the H I and 
stars are coaligned. If the gas and the stars are significantly 
misaligned, then the optically derived inclination angles may 
either overestimate or underestimate the inclination correction 
for the H I line width, depending on the sense of the misalign- 
ment. This then leads to inferred rotational velocities that are 
either too large or too small, with about equal probability. This 
effect must occur at some level, as galaxies with H I disks with 
varying degrees of misalignment are well-known, the most dra- 



matic examples being polar ring systems, in which the gaseous 
and stellar disks are exactly orthogonal to one another. That this 
effect does occur in our sample can be seen in the Tully-Fisher 
diagrams shown in Figure 3: most galaxies cluster around a 
central ridgeline, but a significant number of objects deviate 
both to the low- velocity and high- velocity side of the ridgeline. 
A crucial point to note, however, is that the scatter is not sym- 
metric: there is an excess of low-velocity sources compared to 
high-velocity sources. Moreover, the absolute deviation from 
the Tully-Fisher ridgeline is more extreme for the low-velocity 
objects than the high-velocity objects. Although H I disk mis- 
alignment may account for some of the \ow-v m /ao objects, this 
is not the whole story. 

The asymmetric distribution of the Tully-Fisher outliers 
compels us to conclude that the H I in at least some of these 
objects must reside in a truly dynamically unrelaxed configu- 
ration with respect to the stellar gravitational potential. Three 
possibilities come to mind. First, the H I gas may be strongly 
disturbed as a result tidal interactions, or perhaps distributed in 
tidal tails formed in the aftermath of a galaxy-galaxy merger. 
From optical studies of close pairs and interacting galaxies, 
Barton et al. (2001) and Kannappan & Barton (2004) dis- 
cuss how asymmetric and truncated rotation curves can gener- 
ate strong offsets in the Tully-Fisher relation. While this effect 
might account for some fraction of the outliers, recall that dur- 
ing our sample selection we have purposefully excluded galax- 
ies with close companions and obvious morphological peculiar- 
ities. Thus, this cannot offer a viable solution to the problem at 
hand. 

Second, perhaps the H I owes its kinematic anomaly to a 
large-scale galactic outflow, driven either by strong nuclear ac- 
tivity or a global starburst. We also consider this possibility to 
be untenable. The few galaxies with luminous active galactic 
nuclei or starbursts in our study (e.g., NGC 1068, NGC 3079, 
NGC 405 1) — the most promising candidates for this scenario — 
belong, in fact, to the kinematically normal sample, exactly the 
opposite of what is expected. Furthermore, although we do not 
have uniform data to estimate star formation rates for the en- 
tire sample, the similarity in the broad-band optical and near- 
infrared colors of the two groups strongly suggests that their 
star formation rates are not grossly dissimilar. Lastly, it is worth 
remarking that the amount of neutral gas involved is quite sub- 
stantial. The median H I mass for the low-u m /(7o sample ranges 
from 3 x 10 8 M Q for the SOs to - 1 x 10 9 M Q for the Es, 2 x 10 9 
M Q for Sa spirals, and ~lx 10 10 M & for spirals of type Sb and 
later. Without detailed mapping we do not know precisely how 
much of the gas actually resides in the low-velocity component, 
but for it to markedly affect the integrated profile to the degree 
observed, the fraction of the total mass involved cannot be that 
miniscule. It seems doubtful that nuclear outflows or galactic 
fountains can perturb such a large quantity of gas. In NGC 891 
(Swaters et al. 1997) and NGC 2403 (Fraternali et al. 2002), 
two of the best-studied spirals whose kinematically anomalous 
H I has been interpreted as possibly having a galactic fountain 
origin, the amount of gas affected is ~Tewx 10 8 M Q , roughly 
10% of the total H I mass. Moreover, this material rotates only 
a few tens of km s" 1 slower than the gas in the plane, nowhere 
near as extreme as some of our objects. 

The final and, we believe, most likely possibility is that the 
dynamically unrelaxed H I traces material acquired from an ex- 
ternal source, in the form of a minor merger with a gas-rich 
dwarf galaxy, capture of fallback material from incomplete con- 
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version of gas to stars after a major merger, or cold accretion 
from primordial clouds. Past studies of early-type galaxies have 
long invoked an external origin to account for the unusual prop- 
erties of the neutral hydrogen in these systems, including the 
frequent detection of morphological and kinematic misalign- 
ments between the H I and stars (e.g., Raimond et al. 1981; 
van Gorkom et al. 1986; van Driel & van Woerden 1991; Oost- 
erloo et al. 2002, 2007; Noordermeer et al. 2005; Morganti 
et al. 2006) and the apparent decoupling between the H I and 
stellar mass (Knapp et al. 1985; Wardle & Knapp 1986). In 
this sense, our results for E and SO galaxies could have been 
anticipated. More surprising is that the same phenomenon per- 
sists well into the spiral sequence, to Hubble types as late as 
Scd. The incidence of low-f m /ao objects does decrease toward 
later Hubble types, but this may be an artifact of dilution by the 
intrinsically much larger H I content in these gas-rich systems. 
If, for the sake of argument, external accretion provides a base- 
line supply of 10 8 M & of H I to all sizable galaxies, irrespective 
of Hubble type or environment, this kinematically anomalous 
component would constitute a significantly larger fraction of 
the total H I budget in an E or SO galaxy (~ 20% - 50%) than in 
an Sb or Sc galaxy (~ l%-5%), rendering it much more subtle 
to detect in the latter. 

Notwithstanding the apparent similarity between early-type 
and spiral galaxies, the origin of the kinematically anomalous 
gas may be very different in these two classes of objects. The 
characteristically denser environment around ellipticals and at 
least some SO galaxies sets a natural stage for interaction or 
merger-induced processes to play a more central role than in 
later-type galaxies. Spirals, on the other hand, may have ac- 
cess to another channel of gas supply — "cold accretion." In 
numerical simulations of galaxy formation, accretion of cold 
gas along filamentary structures dominates the growth of lower- 
mass galaxies at high redshifts and in low-density environments 
today (Keres et al. 2005; Maccio et al. 2006). Despite the 
purported cosmological importance of this process, there has 
been very limited observational evidence to support it. A long- 
standing argument that spiral galaxies may accrete appreciable 
amounts of gas directly from the intergalactic medium comes 
from considerations of the high-velocity clouds in the Milky 
Way (Oort 1970; Wakker et al. 1999) and their probable ex- 
tragalactic counterparts (van der Hulst & Sancisi 1988; Thilker 
et al. 2004). Using a blunt but effective tool, namely the v m - 
(Jo diagram in combination with the Tully-Fisher relation, this 
study has highlighted a population of galaxies with kinemati- 
cally anomalous H I that appear to be excellent candidates for 
undergoing cold accretion. Deep aperture synthesis observa- 
tions of isolated spirals are needed to confirm or reject this hy- 
pothesis; to our knowledge, none of the candidates from our 
list has been studied in this way. The selection by environment 
is important in order to distinguish a true cold accretion event 
from other sources of gas infall. For example, the Sc spiral 
NGC 4254, one of the low-u m /tro galaxies in our sample, has 
long been known to possess kinematically distorted H I, but its 
location in the Virgo cluster suggests that its peculiar gas kine- 
matics and distribution most likely arose from a tidal encounter 
with another galaxy or from interactions with the intracluster 
medium (e.g., Phookun & Mundy 1995; Vollmer et al. 2005). 



We reinvestigate the relation between the maximum rotation 
velocity of the disk and the central stellar velocity dispersion of 
the bulge in order to evaluate the claim that these two quantities 
are tightly correlated and independent of galaxy type. Mak- 
ing use of integrated H I line widths to estimate v m and an ex- 
panded database of <jq values augmented by new measurements 
from SDSS, our sample contains 792 galaxies, almost a factor 
of 5 larger than any past study. Contrary to previous reports, we 
find that the v m -ao relation contains significant intrinsic scatter, 
and that the ratio v m /ao varies systematically with Hubble type, 
bulge-to-disk ratio, and light concentration. The density profile 
of a galaxy plays a major role in determining vjao, although 
the large residual scatter suggests that the kinematic structure 
of the galaxy also matters. Extreme late-type spirals that lack 
a clear bulge but contain a central nuclear star cluster deviate 
dramatically from the v, n -<Jo relation. The observed distribu- 
tion of galaxies on the v m -<ro plane serves as an important local 
boundary condition to constrain models of galaxy formation. 

The lack of a tight v m -ao relation removes the principal mo- 
tivation for substituting the halo for the bulge as the galaxy 
component most closely linked with the central black hole. Al- 
though v m is inferior to oo as a predictor of black hole mass, the 
Mbh-v„, relation remains a useful tool in instances when ao is 
too difficult to measure. 

To constrain the intrinsic distribution of v m /ao for any given 
Hubble type, we constructed a /f s -band Tully-Fisher relation 
for the sample using near-infrared photometry from 2MASS. 
The /£j-band Tully-Fisher relation is essentially invariant, for 
morphological types ranging from elliptical galaxies to late- 
type spirals. This exercise revealed an unexpected population 
of outliers characterized by having anomalously low rotation 
velocities for their luminosity, and correspondingly low v m /<Jo 
values. While misaligned H I disks and tidal tails may account 
for some of these low-v m /ao objects, we argue that the major- 
ity of them must have acquired their low-velocity, dynamically 
unrelaxed gas through external capture or cold accretion from 
the intergalactic medium. 
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APPENDIX 

COMPARISON OF LITERATURE DATA 
Stellar Velocity Dispersions 

Among the 293 galaxies from our original sample 1 selected from Hyperleda (see §2), 37 have independent <jq measurements 
given in the SDSS database. Figure 8a compares the velocity dispersion data in common. The quoted error bars from SDSS are 
significantly smaller than those given in Hyperleda, but it is not clear how realistic the SDSS error values are. Although the scatter 
is large, on average the SDSS velocity dispersions tend to be smaller than those given in Hyperleda; for this comparison sample, we 
find ((Jo(SDSS)/(To(Hyperleda)) = 0.89 ± 0.30 (Fig. 8b). Because the Hyperleda velocity dispersion scale has been cross-calibrated 
using multiple measurements of bright galaxies (following the procedure of McElroy 1995), we believe that the Hyperleda values 
are more reliable. Consequently, we have increased the velocity dispersions for the SDSS objects (sample 2) by a factor of 1.12. 

Rotation Velocities 

Reliable integrated H I line widths are available for 93 of the galaxies for which Courteau et al. (2007) compiled rotation velocities 
measured from extended optical rotation curves. A comparison of the two quantities is given in Figure 9a. There is generally 
good agreement between the two sets of measurements, but on average the H I-based velocities are smaller by a factor of 1 .09, 
with a standard deviation of 0.22. Given that our values of v m depend on a specific calibration between H I line width and rotation 
velocity and a specific formalism for correcting the H I line width for turbulent broadening (see §2), this level of agreement is quite 
satisfactory. In any case, only a small fraction of our objects (4%; total of 33 objects from samples 4b and 4c) come from Courteau et 
al.'s compilation. Closer inspection reveals a number of prominent outliers, particularly toward low H I velocities. Figure 9b, which 
plots the ratio of H I to optical velocities as a function of morphological type, illustrates that early-type (r <;0, or E and SO) galaxies 
show the greatest tendency for the H I line widths to underpredict the optical rotation velocities. The implications of this trend are 
further discussed in §4.3. 

We remarked in §2 that Hyperleda's procedure of homogenizing H I line widths does not take into account the possibility of source 
confusion. We illustrate this effect in Figure 10, where the Hyperleda velocities are compared directly with our rederived velocities. 
Although most of the points agree reasonably well for v m ^150-200 km s" 1 , note that at lower velocities the Hyperleda values tend 
to be systematically and significantly (up to ^40%) larger than ours. We attribute this effect to source confusion in the Hyperleda 
average. Since the relative radial velocities of neighboring galaxies are generally larger than the internal velocity of any constituent 
galaxy, source confusion naturally leads to systematically larger line widths. 

Inclination Corrections 

In the course of tracking down the cause of highly discrepant values of v m /<ro, we noticed that some of the inclination angles 
listed in Hyperleda are incorrect. For instance, the nearby, well-resolved late-type spiral NGC 4395 is listed in Hyperleda as having 
b/a = 0.37 or ;' = 90°, whereas the Third Reference Catalogue of Bright Galaxies (RC3; de Vaucouleurs et al. 1991) gives b/a = 0.83, 
which corresponds to ;' = 38°. Simple inspection of the Digital Sky Survey images clearly shows that NGC 4395 is far from edge-on. 
The SO galaxy NGC 3 106 provides another example: Hyperleda lists b/a = 0.30 and i = 90°, while the RC3 gives b/a = 1.0, consistent 
with being nearly face-on, again in agreement with visual examination of the optical image of the galaxy. While most galaxies in 
our sample show less blatant disagreement, we find that there is a general tendency for Hyperleda to underestimate the axial ratio 
(b/a), and hence to overestimate the inclination angle. This is illustrated in Figure 11, where we use the RC3 as a reference. Note 
the systematic difference between the two databases, and that the discrepancy is especially pronounced for barred galaxies. While in 
principle there is no a priori reason to trust the RC3 more than Hyperleda, our systematic inspection of images for our sample (e.g., 
the two cases mentioned above) leads us to believe that the axial ratios given in the former are more reliable than those in the latter. 

It is unclear why the axial ratios listed in Hyperleda are systematically lower. In the case of barred galaxies, which show the 
largest systematic deviations, it is conceivable that the high surface brightness of the bar may have biased the apparent ellipticity of 
the isophotes used to deduce the axial ratio. In some instances, we noted that the final axial ratio adopted in Hyperleda does not, it 
seems, correspond to the B-band data actually tabulated in the database, from which the axial ratio (parameter logr2 5) supposedly 
was derived. Instead, it appears that the adopted axial ratio included in its average size measurements from near-infrared bandpasses. 
In other cases (e.g., PGC 18506), the database does not, in fact, list any fi-band measurements at all, even though logr2 5 strictly 
speaking is said to be derived from B-band data. Since features such as bars tend to be more prominent in redder bandpasses, and 
shallow near-infrared surveys (e.g., DENIS; Paturel et al. 2005) are less effective at picking up the low-surface brightness outer 
regions of the disk component, barred galaxies will generally appear to have a higher ellipticity than if measured from deeper B-band 
photometry. This probably contributes to the effect seen in Figure 1 1 . 

In light of these complications with the Hyperleda axial ratios, we decided to collect our own values from the literature (see Table 
1). Whenever possible, we give preference to the large body of uniform measurements given in the RC3. We use these axial ratios 
to rederive inclination angles, following the same precepts adopted in Hyperleda (Paturel et al. 1997). The inclination angles are 
then used to correct the H I line widths, as well as to estimate internal extinctions for the absolute magnitudes, again following the 
procedures used in Hyperleda. 
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FIG. 8. — Comparison of central velocity dispersions for 37 objects contained in both the Hyperleda sample (sample 1) and the SDSS sample (sample 2). On 
average (o- (SDSS)/cr (Hyperleda)} =0.89 ±0.30. 




100 200 300 400 -4 -2 2 4 6 

v m (H I) [km s" 1 ] T 

FIG. 9. — (a) Comparison of rotational velocities obtained from integrated H I profiles with those measured from resolved optical rotation curves, as compiled for 
93 galaxies in the sample of Courteau et al. (2007). (b) Ratio of rotation velocities derived from H I profile versus resolved optical rotation curves as a function of 
morphological type. Note the strong disagreement for early-type (T < 0) systems. 
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